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Abstract 

We conducted a comprehensive molecular phylogenetic study for a group of 
chameleons from Madagascar (Chamaeleonidae: Calumma nasutum group, com- 
prising seven nominal species) to examine the genetic and species diversity in this 
widespread genus. Based on DNA sequences of the mitochondrial gene (ND2) from 
215 specimens, we reconstructed the phylogeny using a Bayesian approach. Our 
results show deep divergences among several unnamed mitochondrial lineages that 
are difficult to identify morphologically. We evaluated lineage diversification using a 
number of statistical phylogenetic methods (general mixed Yule-coalescent model; 
Speciesldentifier; net p-distances) to objectively delimit lineages that we here con- 
sider as operational taxonomic units (OTUs), and for which the taxonomic status 
remains largely unknown. In addition, we compared molecular and morphological 
differentiation in detail for one particularly diverse clade (the C. boettgeri complex) 
from northern Madagascar. To assess the species boundaries within this group we 
used an integrative taxonomic approach, combining evidence from two indepen- 
dent molecular markers (ND2 and CMOS), together with genital and other external 
morphological characters, and conclude that some of the newly discovered OTUs 
are separate species (confirmed candidate species, CCS), while others should best 
be considered as deep conspecific lineages (DCLs). Our analysis supports a total 
of 33 OTUs, of which seven correspond to described species, suggesting that the 
taxonomy of the C. nasutum group is in need of revision. 



Introduction 

Accurately inventorying and delimiting species is a fun- 
damental prerequisite for numerous fields of biological 
research, such as evolutionary theory, biogeography, and 
conservation planning (e.g., Mayr, 1997). Despite this im- 
portance, biologists are still far from possessing a complete 
species inventory of our planet, which in turn seriously ham- 
pers our ability to conserve biotic resources (UNEP, 2008). 
Even within vertebrates, one of the best-known groups of or- 
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ganisms, there are hundreds if not thousands of species yet to 
be described (e.g., MacKinnon, 2000; Meegaskumbura et al., 
2002; Fouquet et al, 2007; Vieites et al, 2009). The largest 
proportion of undescribed biotic diversity is suspected to 
live in the tropics, and one megadiverse area still containing 
many undescribed species is Madagascar. This, the globe's 
fourth largest island, harbors one of the most spectacular 
but threatened biotas, with many components such as the 
amphibians and reptiles characterized by an extreme degree 
of endemicity. Madagascar's herpetofauna has been more 
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intensively studied in the past two decades compared to most 
tropical areas, yet this scrutiny has paradoxically revealed an 
enormous gap in taxonomic knowledge, and many species 
still await formal description. 

While the largest numbers of undescribed species have 
been identified among amphibians (Kohler et al., 2005; 
Vieites et al, 2009), recent comprehensive taxonomic re- 
views have also revealed that numerous widespread Mala- 
gasy reptile species actually represent species complexes (e.g., 
Raxworthy and Nussbaum, 2006; Gehring et al, 2010, 2011; 
Crottini et al., 2011; Miralles et al, 2011). Although the 
phenomenon of morphologically cryptic species has been 
recognized by taxonomists for nearly 300 years (Bickford 
et al., 2007), it remains largely unexplored among Malagasy 
reptiles. 

The ease of routinely sequencing short DNA fragments 
to identify candidate species (DNA barcoding) has greatly 
facilitated species inventories. However, this represents only 
the first step in integrative taxonomic practice (Padial et al., 
2010), and its uncritical use can lead to poor-quality species 
delimitation because deeply divergent alleles can coexist 
within species due to introgression and incomplete lineage 
sorting, resulting in apparent species paraphyly in mitochon- 
drial DNA (mtDNA; Funk and Omland, 2003). Hence, the 
finding of unexpected genetic diversity in some Malagasy ver- 
tebrates and its translation into species description has led 
to doubts whether current practice represents taxonomic in- 
flation, or is a true representation of species diversity (e.g., 
Isaac et al, 2004; Tattersall, 2007; Weisrock et al, 2010). 
For example, the species diversity of mouse lemurs (genus 
Microcebus) has increased more than sevenfold during the 
last two decades, largely through the analysis of mtDNA se- 
quence data only (Tattersall, 2007; Weisrock et al., 2010). 

Taxonomic progress toward identification and description 
of new species leads to a further challenge in biodiversity stud- 
ies and conservation planning: once a formerly widespread 
species is split into several new species, historical distribution 
records become unreliable with respect to the new taxonomy, 
which subsequently leads to decreased reliability in the bio- 
geographic data available. For example, the splitting of the 
Malagasy chameleon Ccdumma brevicorne into several new 
species (Raxworthy and Nussbaum, 2006) resulted in reliable 
records existing from only a few localities per species (Glaw 
and Vences, 2007). The same situation may hold true for 
the chameleons in the C. nasutum group, which is presumed 
to include six or seven species, widely distributed across the 
rainforests in eastern and northern Madagascar. All of these 
species are easily identified as members of the group by mor- 
phology (most obviously by the soft dermal appendages on 
their snout tips; Fig. 1), and most of them can also be di- 
agnosed to species with some confidence following current 
taxonomy (see Gehring et al., 2011). However, some species 
such as C. nasutum are known to harbor important mor- 



phological variation between and within populations (e.g., 
absence or presence of dorsal crests, shape, length, and col- 
oration of rostral appendages), suggesting a potentially com- 
plex taxonomic situation (e.g., Hillenius, 1959; Brygoo, 1971; 
Glaw and Vences, 2007; Gehring et al, 2011). 

The aim of the present study is a preliminary assessment 
of genetic diversity among these small-sized chameleons. To 
test the hypothesis that previously undiscovered deep genetic 
lineages are present, we (1) produced a phylogeny from a 
fragment of the mitochondrial ND2 gene from 215 samples 
covering the entire distribution of all species, (2) applied 
divergence thresholds using sequence divergence estimates, 
(3) applied a coalescence-based method of species delimita- 
tion (general mixed Yule-coalescent model [GMYC] ) to com- 
pare sequence divergences and subsequently generate candi- 
date species hypotheses, and (4) performed preliminary tests 
of the candidate species hypotheses in one morphologically 
well-delimited complex of these chameleons from northern 
Madagascar (the C. boettgeri complex) based on the inte- 
gration of evidence from mtDNA and nucDNA sequences, 
distribution patterns (i.e., occurrence of syntopy), genital 
morphology, and pholidosis. 

Materials and Methods 
Taxonomy and terminology 

To facilitate the presentation of results we herein use a num- 
ber of predefined terms to refer to units in our phylogenetic 
trees and assemblages of taxa. First, we follow previous au- 
thors, in particular Brygoo (1971), who defined phenetic 
species groups among Malagasy chameleons, and apply the 
term "Calumma nasutum group" to all species in the genus 
Calumma Gray, 1865 with flexible ("nose" like) rostral ap- 
pendages. Species-level taxonomy of these chameleons herein 
follows Gehring et al. (2011), Raxworthy et al. (2008), 
and Andreone et al. (2001), according to which the group 
contains the nominal taxa C. boettgeri (Boulenger, 1888), 
C. fallax (Mocquard, 1900), C. gallus (Gunther, 1877), 
C. guibei (Hillenius, 1959; Fig. 1), C. linotum (Muller, 1924), 
C. nasutum (Dumeril and Bibron, 1836), and C. vohibola 
(Gehring et al., 2011). Within the C. nasutum group, three 
species from northern Madagascar (C. boettgeri, C. guibei, 
C. linotum) Are easily and immediately recognized by having 
well-defined occipital lobes, although it is uncertain whether 
these species form a monophyletic group. We refer to these 
three species herein as the C. boettgeri complex. 

The species delimitation approaches used herein are based 
on a phylogenetic tree derived from an mtDNA gene, and 
species hypotheses were produced by drawing species limits 
at certain nodes of the tree. We refer to genealogical units 
at or below the species level as lineages. Those lineages sup- 
ported by the species delimitation approaches (see below) 
are here termed operational taxonomic units (OTUs). For a 
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Figure 1. Male specimen of Calumma guibei from the Tsaratanana 
Massif, Madagascar, in life. Note the rostral appendage that is flexi- 
ble and dermal (no underlying bony structure) in all species of the C. 
nasutum group, and the small occipital lobes which are characteristic for 
the species of the C. boettgeri complex. 



subset of species, additional evidence from morphology or 
nuclear genes was available; following Vieites et al. (2009) 
we use the terms confirmed candidate species (CCS) for lin- 
eages supported by at least two independent lines of evidence, 
and deep conspecific lineage (DCL) if the available evidence 
rejects evolutionary independence and thus species status. 

Some of the lineages identified are morphologically sim- 
ilar and are grouped phylogenetically to form a number of 
well-supported major clades. Furthermore we use the term 
"cluster" in some cases to refer to the results of some of the 
species delimitation approaches which are based on cluster- 
ing sequences on the basis of their similarity (see below). 

Sample collection and morphology 

Specimens were collected in the field largely at night using 
torches and headlamps to detect roosting chameleons in the 
vegetation. Voucher specimens were anesthetized and euth- 
anized, fixed with 90% ethanol and stored in 70% ethanol. 
Muscle tissue samples for molecular analyses were taken from 
the tongue, hind limb, or tail and preserved in pure ethanol. 
Individuals that were not collected as vouchers were sampled 
and released. All tissue samples were stored in 99% ethanol. 
ZCMV, FGZC, DRV, MPFC, ACZC, ZCSH, and PSG refer to 
field numbers of M. Vences, F. Glaw, D. Vieites, M. Pabijan, 
A. Crottini, S. Hauswaldt, and P.-S. Gehring, respec- 
tively. Representative voucher specimens were collected 
and deposited in the herpetological collections of the 
Museo Nacional de Ciencias Naturales (MNCN), Zo- 
ologisches Forschungsmuseum Alexander Koenig, Bonn 
(ZFMK), Zodlogisch Museum Amsterdam (ZMA), Zool- 
ogische Staatssammlung Miinchen (ZSM), and the Uni- 
versite d' Antananarivo, Departement de Biologie Animale 



(UADBA). Geographic coordinates were recorded with GPS 
receivers. (Table SI provides detailed information on GPS 
coordinates of sampling locations.) A distribution map in- 
cluding an elevation layer of sample localities in Madagas- 
car was obtained using the program DIVA GIS (Version 
7.3.0.1; http://diva-gis.org). The taxonomy used herein fol- 
lows Raxworthy et al. (2008) for the definition of C. boettgeri, 
C. guibei, and C. linotum, and Glaw and Vences (2007) and 
Gehring et al. (2011) for C. fallax, C. nasutum, and C. vohi- 
bola. Samples of C. oshaughnessyi (FGZC 4577) and Brookesia 
ramanantsoai (FGZC 4557) were used as outgroups in the 
molecular analyses. 

For study of hemipenes, these structures were everted di- 
rectly after euthanasia, by first injecting water and subse- 
quently fixative into the organs. Hemipenes were examined 
under a stereo microscope and hand drawn. Variation and 
possible artifacts caused especially by incomplete eversion 
of organs or detailed structures were assessed by examin- 
ing multiple specimens per lineage in the C. boettgeri com- 
plex. Terminology of hemipenial structures follows Klaver 
and Bohme (1986) and Bohme (1988). 

Laboratory protocols 

Total genomic DNA was extracted from the tissue sam- 
ples using proteinase K digestion (10 mg/mL concen- 
tration) followed by a salt extraction protocol (Bruford 
et al., 1992). A fragment of the mitochondrial gene 
for NADH dehydrogenase subunit 2 (ND2) was PCR- 
amplified using the primers ND2F17 (5'-TGACAAAAAAT 
TGCNCC-3') (Macey et al, 2000) and ALAR2 (5' -AAA 
ATRTCTGRGTTGCATTCAG-3') (Macey et al, 1997), and 
a fragment of the nuclear DNA (nucDNA) from the 
oocyte maturation factor gene (CMOS) using the primers 
C08 (5'-GCTTGGTGTTCAATAGACTGG-3') and C09 (5'- 
TTTGGGAGCATCCAAAGTCTC-3' ) (Han et al, 2004), fol- 
lowing standard protocols. After purification with ExoSAP- 
IT (USB), the fragments were sequenced with an auto- 
mated DNA sequencer (ABI 3130 XL Applied Biosystems, 
Foster City, USA). For the ND2 dataset only the light 
strand was sequenced, while the CMOS gene fragment 
was sequenced in both directions to verify possible het- 
erozygotes. Quality of sequence chromatograms was verified 
and sequences were aligned manually by amino-acid transla- 
tions with the software CodonCode Aligner (CodonCode 
Corporation, Dedham, USA). All newly determined se- 
quences were submitted to GenBank (accession numbers 
JQ733989-JQ734169 and JQ614354-JQ614415; see also Ta- 
ble SI). 

Phylogenetic analysis 

Three sets of DNA sequences were used for analyses: ( 1 ) a 
586 bp fragment of ND2 from all available samples (215 
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specimens; 103 unique sequences), (2) a subset of these 
same ND2 sequences trimmed to exactly the same length 
in all taxa (389 bp; no missing data) from only chameleons 
with flexible rostral appendages as well as occipital lobes 
(the C. boettgeri complex; N = 82), and (3) a 500 bp frag- 
ment of CMOS for most samples of the C. boettgeri com- 
plex (N = 63). The major focus of the present paper is on 
species delimitation, not on detailed reconstruction of the 
phylogeny of the C. nasutum group, for which the short 
DNA fragments used herein would be clearly insufficient. 
We analyzed three separate datasets using multiple com- 
plementary methods to examine the hypothesis that deep, 
potentially species-level divergences exist within the species 
complex: 

(1) Phylogenetic analysis of the full ND2 datasets with each 
unique sequence represented only once (N = 103) was car- 
ried out using Bayesian inference in BEAST (Drummond 
and Rambaut, 2007), in order to obtain an ultrametric 
tree for automated species delimitation analysis (see sec- 
tion "Species Delimitation" below for details). In order to 
graphically represent the phylogenetic placement of the se- 
quences of all specimens (N = 215) we also ran a Bayesian 
inference analysis with MrBayes version 3.1.2 (Ronquist and 
Huelsenbeck, 2003) under the same settings as for the 
C. boettgeri complex data (given in the next paragraph); 
this tree is only represented in the Supporting information 
(Fig. SI) and was not used further in the study. 

(2) Analysis of the subset of ND2 sequences of only the 
C. boettgeri complex was carried out with MrBayes un- 
der substitution models determined with MrModeltest 2.3 
(Nylander, 2004) with two independent runs of 20 million 
generations, started on random trees and each with four in- 
crementally heated Markov chains (using default heating val- 
ues), sampling the Markov chains at intervals of 1000 gener- 
ations. Convergence of the Markov chains was checked using 
Tracer vl.5 (Rambaut and Drummond, 2009) and mixing of 
chains was assessed with AWTY (Nylander et al., 2008). The 
first 25% of the trees were conservatively discarded and the 
trees of both runs retained post burn-in were used to generate 
a 50% majority rule consensus tree. 

(3) CMOS sequence was used as independent genetic marker 
to confirm the absence or reduction of gene flow among ge- 
nealogical lineages initially identified on the basis of the ND2 
analysis. Because haplotype sharing among lineages can also 
reflect incomplete lineage sorting, this kind of evidence is 
asymmetric: haplotype sharing in independent markers can- 
not reject reduction of gene flow but two groups of specimens 
not sharing haplotypes of multiple independent markers are 
likely to be on independent evolutionary trajectories without 
gene flow (the genealogic concordance method of phyloge- 
netic species recognition, proposed by Avise and Ball, 1990). 
Because we found only minor variation in CMOS, we used 
a haplotype network approach for this gene. Sequences were 



phased into distinct haplotypes using the Phase algorithm 
implemented in the software DNASp version 5.10.3 (Librado 
and Rozas, 2009). A haplotype network was then constructed 
from the nuclear haplotypes using TCS (Clement et al, 2000). 
Ambiguities in the network were solved manually using fre- 
quency, topological, and geographical criteria (e.g., Posada 
andCrandall, 2001). 

Species delimitation 

To examine the genetic diversity in this group, three ap- 
proaches were used on the basis of the full set of ND2 se- 
quences (506 bp, 215 specimens; 103 unique sequences). 
Of these the approaches 1 and 2 (see below) searched for 
the breaks in overlap between intra- and interspecific vari- 
ation, whereas approach 3 was a coalescent-based method 
that attempts to identify boundaries between species-level 
and population-level processes from tree shape and branch 
lengths. 

In approach 1, the topology obtained in the Bayesian anal- 
ysis was used as a framework for manually grouping taxa at 
different levels (excluding duplicate haplotypes, remaining 
sequences for analysis N = 103) by searching for the thresh- 
old at which the intraclade p-distances in the ND2 sequences 
were lower than the interclade p-distances. We performed 
this approach using net p-distances, that is, the average p- 
distances between lineages subtracting the differences within 
each of the lineages. Taxa were grouped into clades start- 
ing at the base of the tree, and moving outward toward the 
tips until supported nodes (PP of 0.95 or higher) were en- 
countered. The net p-distances among and within these deep 
lineages were then calculated. The taxa were then regrouped 
into different combinations of lineages by moving outward 
to the next set of supported nodes and the net p-distances 
estimated again. The procedure was repeated a final time, 
moving outward to the next set of supported nodes, and 
p-distances again estimated for these shallow lineages. 

In approach 2, the Speciesldentifier "Cluster" algorithm 
in Taxon DNA v 1.7 (Meier et al, 2006) was used to au- 
tomatically cluster taxa according to pairwise distances for 
every sequence within each cluster. Incremental threshold 
values of 2-10% were used, whereby the maximum pairwise 
distance within each cluster (possible species-level lineage) 
should not exceed a given threshold. 

Finally, in approach 3, clustering was examined using a 
GMYC in which the boundaries between species-level and 
population-level processes are identified (Pons et al., 2006; 
Fontaneto et al., 2007). In this method, the transition in 
branching rates from lineage diversification processes (Yule 
model) to population-level diversification (coalescent model) 
are identified. The GMYC model can be employed with a sin- 
gle transition (Pons et al., 2006) or a multiple transition that 
allows the depth of the coalescent-speciation transition to 
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vary across the phylogenetic tree (Monaghan et al, 2009). 
The observed branching rate is then compared against a 
null model that assumes no shift in branching rate using 
a log-likelihood ratio test (LRT). Failure to reject the null 
model implies no shift in the branching rate (i.e., a single 
species is present). An ultrametric tree was created in BEAST 
(Drummond and Rambaut, 2007), using a Yule model for the 
tree prior, a strict molecular clock (fixed to a rate of 1.0), a 
random starting tree, and the GTR+I+G model of sequence 
evolution. Tracer v. 1.5 (Rambaut and Drummond, 2009) 
was used to verify that the effective sample size for all es- 
timated parameters exceeded 200. The ultrametric tree was 
obtained using TreeAnnotator in BEAST (Drummond and 
Rambaut, 2007) specifying a 10% burn-in, the maximum 
credibility clade tree, and saving the target node heights. To 
ensure that the choice of priors in BEAST did not affect 
the result, the ultrametric tree was also obtained using ( 1 ) a 
Yule model with a lognormal relaxed clock, (2) a coalescent 
model with a strict clock, and (3) a coalescent model with a 
log-normal relaxed clock, following Monaghan et al. (2009). 
The four trees obtained were fit to both single-transition 
and multiple-transition GMYC models using the Splits and 
Ape packages in R. The fit of the single and multiple GMYC 
models were compared with a / 2 -test in R using the Splits 
package (Monaghan et al. 2009). 

Focusing mainly on the C. boettgeri complex, we further- 
more followed an integrative taxonomy approach that uti- 
lizes the best available evidence from different methods and 
approaches (Padial et al, 2010). We took (1) syntopic oc- 
currence without genetic admixture and with maintenance 
of morphological differences, (2) congruent signals from 
largely independent character sets (morphology, mtDNA, 
nucDNA), and (3) distinct and constant differences in char- 
acters mediating reproductive isolation (i.e., genital morpho- 
logy) as indications of evolutionary independence and thus 
species status of lineages. 

Results 

ND2 phylogeny and sequence diversity 

Altogether, 103 unique ND2 sequences were obtained from 
215 Calumma individuals from more than 60 different loca- 
tions (Fig. SI and Table SI), and after exclusion of duplicate 
haplo types, 103 sequences were used to create the ultrametric 
tree in BEAST (Fig. 2; a detailed version of this tree is given 
in Fig. S2). 

Based on net p-distances for the entire ND2 dataset 
(pDist-approach 1; Fig. SI), three sets of possible species- 
level lineages were identified: (1)11 lineages with divergences 
>8.2% between and <8. 5% within them; (2) 27 lineages with 
divergences >6% between and <4.5% within them, and 
(3) 32 lineages with divergences >6% between and <3% 
within them (Fig. 2; Table S2). Because the first set of diver- 



gences contained an overlap in values for the between and 
within p-distances, this set was considered the least likely to 
accurately reflect the number of independent evolutionary 
lineages. Both the second and third sets had gaps in values 
for the between and within ^-distances, and in both cases the 
divergence among lineages was at least 6%. However, for the 
third set, the within-lineage divergences were smaller (<3%) 
than those within the second set ( <4.5%), resulting in alarger 
gap for the values within and between lineages. 

Using Speciesldentifier (SpeciesID — approach 2), the 
number of species-level lineages identified depended on 
the threshold used, with the largest threshold producing the 
fewest lineages (Table S3). Given that approach 1 shows that 
between-lineage p-distances are maximized when the within- 
lineage value is no more than 4.5%, a threshold of 4% or 5% 
should produce similar results for approach 2. For these data, 
the 5% threshold produced 34 lineages, a value similar to that 
of approach 1 . 

The number of species-level lineages defined by the GMYC 
model (approach 3) relied somewhat on the priors used to 
construct the ultrametric tree, ranging from 44 to 53 OTUs 
(Table S4). There was no significant difference between the 
single- and multiple-threshold approaches for any tree mod- 
els. The single threshold GMYC model utilizing a tree based 
on a coalescent model and a strict clock produced the least 
number of clusters (23) but not the fewest singletons (24). 
Overall, the Yule model with the strict clock as well as the 
coalescent model with the relaxed clock produced fewer to- 
tal lineages (44), comprised of 26 clusters and 18 singletons. 
For further analysis, we applied a conservative approach and 
identified those lineages that were supported by at least two 
out of three methods (p-distance criteria, Speciesldentifier, 
and GMYC). This yielded a total of 33 OTUs, which we as- 
signed to species complexes based on morphological char- 
acters considered to be diagnostic among species as pre- 
viously defined (mainly presence vs. absence of occipital 
lobes, and size and shape of rostral appendage) and place- 
ment of topotypical populations of some species in the tree 
(Fig. 2). 

Mitochondrial and nuclear sequence 
diversity in the C. boettgeri complex 

Bayesian inference was carried out on ND2 sequences from 
82 individuals of chameleons of the C. nasutum group with 
distinct occipital lobes, assigned to C. boettgeri, C. guibei, and 
C. linotum. The alignment of 389 nucleotides had 211 variable 
sites of which 153 were parsimony informative. The analysis 
(Fig. 3a; duplicate sequences not removed) confirmed within 
this complex the same major lineages as the analysis of the 
full dataset of 103 unique sequences in Figure 2. Specimens 
from Montagne d'Ambre in northern Madagascar formed 
clades corresponding to higher and lower elevations on this 
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Figure 2. An ultrametric tree for the ND2 gene fragment of Malagasy chameleons of the Calumma nasutum group, using Bayesian inference analysis 
with BEAST, a Yule model, and a strict clock. Calumma oshaughnessyi and Brookesia ramanantsoai were used as outgroup taxa (not shown). In the 
tree on the left Bayesian posterior probability values above a threshold of 0.9 are indicated by broadened lines; strongly supported major clades are 
indicated in colors and numbered with capital letters. Different symbols characterize lineages that were flagged as candidate species by two of three 
different methods of species delineation. Nomina currently considered as valid are tentatively assigned to these lineages based mainly on their type 
localities. Locality information for each lineage is given in the map of Madagascar. In the ultrametric tree on the right the lineages are indicated based 
on three different clustering methods (net p-distances: black dots, Speciesldentifier: gray bars; and GMYC clusters: branches in red). Black stars denote 
groups for which all three methods were in agreement, white stars denote groups for which net p-distances and Speciesldentifier were in agreement, 
and gray stars indicate groups for which GMYC and Speciesldentifier were in agreement. 
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N= 3 7 3 1 5 2 




Figure 3. Differentiation among clades in Malagasy chameleons of the Calumma boettgeri complex (species of the C. nasutum group with occipital 
lobes: C. boettgeri, C. guibei, C. linotum). (a) Majority-rule consensus tree from a Bayesian inference analysis based on a fragment of the ND2 gene, 
with exemplary drawings of the hemipenis (asulcal view left, sulcal view right) for major lineages. For convenience, roman numbers were assigned 
to each lineage. Asterisks denote support from Bayesian posterior probabilities: *PP 0.95-0.98; **PP 0.99-1.0. Abbreviations CCS, UCS, and DCL as 
in Figure 1 . (b) Haplotype network based on a fragment of the nuclear CMOS gene (400 bp; 126 haplotype sequences from 63 specimens). Sizes of 
circles on the left indicate haplotype frequencies, (c) Differentiation of lineages in the number of large scales on the nasal appendage of males (dots, 
boxes, and whiskers represent mean, standard error, and minimum-maximum values, respectively); note the highly divergent value in lineage II of 
C. guibei. (d) Distribution map of localities sampled for this study (which cover the entire known distribution area of these species). Colors and symbols 
in map as in Figure 2; bo, gu, and li refer to C. boettgeri, C. guibei, and C. linotum, with lineage numbers following the abbreviation. 
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mountain, suggesting a possible adaptive diversification. The 
haplotype network derived from CMOS sequences of 63 
individual chameleons (duplicate sequences not removed; 
126 haplotype sequences after phasing) showed a clear sep- 
aration of the three species. Samples assigned to clade D 
were separated by a minimum of five mutational steps from 
those assigned to the clades E and F, while samples of those 
two clades were separated from each other by a minimum of 
three mutational steps (Fig. 3b). Nuclear haplotype sharing 
was observed between lineages I and II of clade E, and be- 
tween lineages II and III of clade D, but not between lineages 
I and II/III of clade D (roman numbers as used in Fig. 3). The 
haplotype found in the single specimen of lineage II of clade 
F was separated by 2-3 mutational steps from the haplotypes 
of lineage I of that clade (C. guibei). 

Differentiation in genital and external 
morphology in the C. boettgeri complex 

Hemipenes were examined from a total of 16 adult males of 
the C. boettgeri complex. The hemipenes of the three clades 
representing the complex (clades D, E, and F; Fig. 3a) could 
be diagnosed by various differences, such as the apical pap- 
illary field on the sulcal side present only in clade A and the 
apical sulcal lobes present only in clade F (Table 1). Several 
important differences were also found among lineages within 
clades. The single available male of lineage II in clade F had 
much smaller and shallower calyces and larger lateral rotulae 
as compared to the various males examined from lineage I in 
clade F (C. guibei). In clade E, the development of the sulcal 
lips appears to be different between the two main lineages. 
And in clade D, the papillary field might be more strongly de- 
veloped in lineage I (C. boettgeri) than in lineage II, although 
this possible difference requires more in-depth study. 

The shape and number of scales on the rostral appendage 
in males obviously differs among various species in the 
C. nasutum group (e.g., very long and pointed in C. gallus 
and related species; Fig. 2). Within the C. boettgeri complex 
we noted that the appendage is made of numerous large 
scales, with very tiny intercalated scales that are difficult to 
count. In lateral view, 10-47 small intercalated scales were 
counted in all individuals of clades E and F (N = 4 and 7), 
but they were missing in most individuals of the clade D 
(JV = 10; 0-6 scales in eight individuals; 18-24 scales in two 
individuals). The number of large scales laterally visible on 
the appendages was 43-70 except in the single individual of 
lineage II of clade F. This specimen, characterized by a par- 
ticularly long appendage, had 103 large scales laterally visible 
on its appendage. Small sample sizes made statistical analyses 
of these differences impossible. Pending study of additional 
material, more morphometric and meristic data and detailed 
hemipenial descriptions will be published in forthcoming 
taxonomic revisions of these chameleons. 
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Discussion 

Mitochondrial lineages and their 
distribution pattern 

Our analyses based on the ND2 dataset revealed 33 mito- 
chondrial lineages within the C. nasutum group that were 
supported by at least two of the three species delimitation 
approaches (Fig. 2). These congruent results from multi- 
ple species delimitation methods suggest that current taxon- 
omy may greatly underestimate species diversity within this 
group. 

At deeper levels in the tree, we can further distinguish 
1 1 well-supported clades with net p-distances in excess of 
8.5% between them. These 11 deeper groupings are largely 
congruent with obvious morphological similarities within 
these clades and differences between them. Morphologically 
they roughly agree with current species-level classification: 
clade A corresponds morphologically to C. gallus, clade C to 
C. vohibola, clade D to C. boettgeri, clade E to C. linotum, 
clade F to C. guibei, clades H and I to C. fallax (recovered 
as monophyletic group but without statistical support), and 
clades B, G, J, and K to C. nasutum (Figs. 2 and SI). For 
example, similarity of all specimens in clade A to C. gallus is 
based on their relatively long and pointed rostral appendages. 
Besides the nominal species C. gallus from central Madagas- 
car, five OTUs exist within this clade, all of which inhabit the 
lowland east coast below 900 m (Fig. 2). Males of several of 
the six OTUs within clade A show differences in their ros- 
tral appendage (such as spiny vs. smooth scales). The length, 
shape, and coloration of rostral appendages are thought to 
be important visual signals most probably influenced by sex- 
ual selection via mate choice by females (e.g., Parcher, 1974; 
Gehring et al., 2011), although it is uncertain whether the 
characters showing differences between species are relevant 
to intraspecific communication. 

Clades H and I contain specimens morphologically similar 
to G fallax, although this species is morphologically poorly 
defined in comparison with C. nasutum. These two clades 
form a monophyletic group and are mainly distributed at 
middle and higher elevations. Clade H contains six deep mi- 
tochondrial lineages (C. fallax and five OTUs), all of which 
were sampled from mid-elevations (above 800 m a.s.L; many 
from above 1200 m a.s.L). Clade I includes three highly diver- 
gent OTUs from the Tsaratanana massif in northern Mada- 
gascar. Our analyses found four major clades with morpho- 
logical resemblance to C. nasutum, with a distribution mainly 
at lower elevations (below 800 m a.s.L). Clades D-F (the 
C. boettgeri complex) appear to be restricted to northern 
Madagascar, with a clear elevational zonation: Clade D at low 
to mid-elevations (0 to ca. 1300 m a.s.L), clade E at inter- 
mediate elevations (1300-1550 m a.s.L), and clade F at the 
highest elevations recorded for any species of the C. nasutum 
group (1550-2000 m a.s.L). 
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Table 1. Summary of hemipenial characters observed within lineages of the C. nasutum group characterized by distinct occipital lobes (C. boettgeri, 
C. guibei, C. linotum). Hemipenial structures that are unique within one species clade are marked with one asterisk. Reliable diagnostic characters 
distinguishing one lineage from the other lineage within the same clade are in italics and marked with two asterisks. 



Lineage 


Calyces 


Sulcal lips 


Rotulae 


Elongated papillae 


Apical papillary 
field 


Apical sulcal lobes 


C. aff. boettgeri 
clade D, OTU 8 

C. boettgeri clade 
D 


Large, deep 
(sulcal, asulcal) 

Large, deep, 
(sulcal, asulcal) 


Well developed 
Well developed 


Large sulcal, small 
asulcal pair 

Large sulcal, small 
asulcal pair 


Laterally between 

sulcal and 

asulcal rotulae 
Laterally between 

sulcal and 

asulcal rotulae 


•Medially between 
rotulae 

•Medially between 
rotulae, well 
developed 


Absent 
Absent 


C. guibei clade F 

C. aff. guibei clade 
F, OTU 1 2 


"Large, deep, 

angular (sulcal, 

asulcal) 
"Small, rather 

elongated, 

shallow (sulcal) 


Well developed 
Thin 


** Two pairs of 
small rotulae 

"Small distal, 
relatively larger 
lateral rotulae 


Laterosulcal 
position 

Between distal 
and lateral 
rotulae 


Absent 
Absent 


•Present 
•Present 


C. linotum clade I 

C. linotum clade I, 
OTU 10 


Large, shallow 
(sulcal) 

Large, shallow 
(sulcal) 


** Well developed 
"Relatively thin 


•Large sulcal, small 
distal pair 

•Large sulcal, small 
distal pair 


Between sulcal 

and apical 

rotulae 
Between sulcal 

and apical 

rotulae 


Absent 
Absent 


Absent 
Absent 



Candidate species and prospects for future 
research 

Two of the species delimitation approaches used herein rely 
on sequence thresholds, and we chose a threshold of 6% to 
define OTUs. This value is not unrealistically low; in other 
groups of chameleons (i.e., Bradypodion, Kinyongia) species 
are characterized by divergences ranging from 5% to 25% for 
ND2 (Tolley et al. 2004, 2006, 20 1 1 ) . However, it is clear that 
using a threshold near the lower limit of divergence values 
observed between well-established chameleon species implies 
a risk of overestimating true species numbers. If this evidence 
were to be uncritically applied to the C. nasutum group, and 
all OTUs accepted as candidate species, our results would 
suggest a more than fourfold increase in species numbers, 
from 7 to 33, which probably is an overestimate rather than 
underestimate of the true number of species in the group. 

More detailed (although still preliminary) data from gen- 
ital morphology, pholidosis, and nucDNA are available for 
clades D-F forming the C. boettgeri complex. We suggest 
that clades D-F contain three valid nominal species, two 
unnamed independently evolving units representing species 
(CCS), but also three examples of lineages that are best inter- 
preted as intraspecific variation (DCL). At present, we have 
no morphological or genital morphological data available for 
other clades in the C. nasutum group. However, nuclear al- 
leles are not shared among several lineages morphologically 
similar to C. nasutum (Gehring et al., 2011), and our ob- 



servations (Table 2) indicate various instances of syntopic 
occurrence of lineages with apparent maintenance of mor- 
phological distinctness. This indicates that numerous other 
mtDNA lineages represent undescribed species, even if the 
total species diversity in the group will probably stabilize at a 
number lower than 33. 

Calumma nasutum was assumed to be a widespread species 
in Madagascar's eastern rainforest band (Brygoo, 1971). Our 
data show a high genetic subdivision of this and all other 
species in the C. nasutum group, and if a large proportion 
of the mtDNA lineages indeed represent species, then almost 
no species of the group would have a range spanning over 
more than ca. 200 km in a north-south direction, similar 
to other groups of microendemic organisms in Madagascar 
(e.g., Wollenberg et al, 2008). 

Clearly, our results indicate that a systematic revision and 
the description of various new species within the C. nasutum 
group are necessary. Such a revision will require intensive 
work because various important questions remain unclear: 
(1) monophyly of the C. nasutum group cannot be con- 
firmed without the addition of taxa from other closely re- 
lated species in a phylogenetic context, and including longer 
DNA sequences from additional genes. (2) Testing with a 
comprehensive nuclear dataset the monophyly of the C. 
boettgeri complex. These species (C. boettgeri, C. guibei, and 
C. linotum) are the only C. nasutum group chameleons 
with distinct occipital lobes, are otherwise morphologically 
similar to each other, and are restricted to a small area 
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Table 2. Table of localities in which different lineages of the Calumma nasutum group were found in direct syntopy (less than 1 km from each other). 



Locality 


Direct syntopy of 


Discreteness of lineages evidenced by 


Sahafina 


C. gallus (clade A) 

C. aff. nasutum (clade G; OTU 16) 


(1) Mitochondrial differentiation; (2) constant differences in external 
morphology of males and females 


Anosibe An'Ala 


C. aff. nasutum (clade B, OTU 7) 

C. nasutum (clade K) 

C. aff. fallax (clade H, OTU 20) 


Mitochondrial differentiation 


Tarzanville 


C. aff. nasutum (clade B, OTU 7) 

C. nasutum (clade K) 

C. aff. gallus (clade A, OTU 3) 


(1) Mitochondrial differentiation; (2) constant differences in external 
morphology of males and females (only between C. aff. gallus and the 
two C. aff. nasutum) 


Ambatofotsy 


C. aff. nasutum (clade B, OTU 7) 
C. aff. gallus (clade A, OTU 3) 


(1) Mitochondrial differentiation; (2) constant differences in external 
morphology of males and females 


Makira forest, Campsite I 
(Angozongahy) 


C. aff. nasutum (clade A, OTU 6) 
C. aff. nasutum (clade G, OTU 14) 


(1) Mitochondrial differentiation; (2) differences in external morphology of 
males 


Tsaratanana, Campsite I 
(Antevialambazaha) 


C. guibei (clade F) 

C. aff. linotum (clade E, OTU 1 0) 


Morphological differences of males (no linotum sequence from this site 
available, but specimens observed in the wild) 


Bemanevika forest, Campsite II 
(near village) 


C. guibei (clade F) 
C. linotum (clade E) 


(1) Mitochondrial differentiation; (2) differentiation in nuclear gene; 
(3) morphological differences of males 


Bemanevika Campsite I 


C. linotum (clade E) 

C. aff. //nofum (clade E, OTU 1 0) 


Only mitochondrial differentiation; no differences in morphology, genital 
morphology, or nuclear DNA 



of northern Madagascar (Fig. 3), but the ND2 phylogeny 
(Fig. 2) does not support their monophyly. (3) The iden- 
tity of the available nomina, including several junior syn- 
onyms, needs to be clarified. For instance, we here fol- 
lowed the definition of Raxworthy et al. (2008) who used 
the name C. linotum for populations of chameleons from 
mid-elevations in the Tsaratanana massif, and C. guibei for 
populations at higher elevations in the same regions, but 
since the type locality of C. linotum is unspecified and the 
holotype of C. guibei is a juvenile specimen, this definition 
requires confirmation. Similarly, the identity of C. gallus and 
C. nasutum each require further study (Gehring et al., 201 1). 

In summary, the present study strongly suggests that the 
diversity in the Malagasy chameleon fauna is far higher than 
previously realized. Clearly this understudied fauna merits 
consideration in conservation efforts, since some of these 
newly discovered lineages are probably restricted to small 
geographic areas and face serious danger of extinction due to 
the rapid ongoing destruction of primary habitats. 
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